Nuclear thermal to electric power conversion carries the promise of longer duration 
missions and higher science data transmission rates back to Earth for a range of missions, 
including both Mars rovers and deep space missions. A free-piston Stirling convertor is a 
candidate technology that is considered an e_cient and reliable power conversion device for 
such purposes. While already very e_cient, it is believed that better Stirling engines can 
be developed if the losses inherent in current designs could be better understood. However, 
they are di_cult to instrument and so e_orts are underway to simulate a complete Stirling 
engine numerically. This has only recently been attempted and a review of the methods 
leading up to and including such computational analysis is presented. And _nally it is 
proposed that the quality and depth of Stirling loss understanding may be improved by 
utilizing the higher _delity and e_ciency of recently developed numerical methods. One 
such method, the Ultra HI-FI technique is presented in detail. 
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Stirling Convertor CFD 

System Overview 






Contours of Velocity Magnitude (rn/s) (Time=0.000Qe+00) Mar 22, 2004 

FLUENT 6.1 (axl, dp, segregated, dynamesh, ske, unsteady) 










Schematic with Springs and Dampers 



Effect of Heating & 



/ Bearing System 

Moving Mesh Loss 



Cooling / Artificial Transient Heat Loss 





Flow Losses (or Artificial CFD Losses) 


Gas Spring 



Piston/Displacer Coupling via Flow Losses 



Load 


Working Fluid (bounce space and working space) 



Martini Nomenclature (West) 

(See also Organ and Urieli) 

1. First Order (ideal loss-free analysis, use correction factor) 

- Schmidt or Cooke-Yarborough 

2. Second Order (decoupled losses, use adiabatic analysis) 

- Berchowitz, D.M. (includes dynamics) 

• Subtract heat losses 

- Appendix gap (shuttle and pumping) 

- Regenerator imperfections 

• Subtract power losses 

- Pressure drop (non-uniform flow across regenerators) 

- Transient heat transfer 

3. Third Order (control volumes or nodes, direct solve 1 D equations ) 

- Finkelstein (first), Urieli, Berchowitz (includes dynamics) 

- Implicit Space and Time (GLIMPS, SAGE, finite difference, 1 D) 

- Linearized harmonic analysis (HFAST, finite volume, 1 D) 
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Fourth Order? Full CFD Analysis (2D or 3D) 




Flow distortion (oscillating flow distorts velocity and temperature 
profiles) in the regenerator has significant effects on performance 
and must be modeled with at least two dimensions (Berchowitz). 

• Manifest (2D, D. Gedeon) 

- Beam and Warming 

- Uses thermal non-equilibrium porous media model 

• Fluent (2D and 3D) 

- SIMPLER and PISO 

- Uses thermal equilibrium porous media models 
» Modified CAST (Tew, Ibrahim, Peric, et. al.) 

- SIMPLE 

• CFD-ACE (2D and 3D) 

- SIMPLE & SIMPLER 

- Uses thermal equilibirium porous media models 

• Star-HPC, CFX 

- SIMPLE & PISO 


Whole Engine Simulation 



• Mahkamov, et. al, 
University of Durham 

• Ibrahim, et. al., Cleveland 
State University 

• Wilson, et. al., NASA GRC 

• Recent proprietary 
commercial attempts 




Contours of Velocity Magnitude (m/s) (Tlme=Q.00Q0e+GQ) Mar 22, 2004 

FLUENT 6.1 (axl, dp, segregated, dynamesh, ske, unsteady) 


ALL USE 2 nd ORDER SIMPLE/PISO RANS APPROACH 


TOO SLOW! 
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Too Much Numerical Error 
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Turbulence Modeling for Oscillating Flows 



Incorrect Pressure Drop and Heat Transfer 


• Wavelet Analysis of Multiscales 

• DNS with Wall Layer Approximations 

• Detached Eddy Simulation (RANS/LES) 

• Cascade Technologies, V2F for low Reynolds number 
heat transfer 

• Ultra-wave Kolmogorov Scale Resolution Techniques 


Full DNS Jet Simulations with 
Re>1 00,000 are now possible 


Parallel Solution 
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•Assume 20 micron wire diameter, 25mm by 5mm by 5mm 
regenerator region to simulate 

•Might have (at most) 312000 wires in 2D or 8 million in 3D 
depending on porosity 

•Regularized screen meshes may enable more detailed 
analysis 




120 Million Cells, 128 processors, 24 hours 


Park, J. 


STAR-HPC 



Curved Spectral Volumes (Z.J. Wang) 
Entropy Production and Density Contour 
M=0.3, Inviscid Flow 



2 nd Order Piecewise Linear 2 nd Order “High Fidelity Grid” 3 rd Order Piecewise Quadratic 
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Curved Spectral Volumes (Z.J. Wang) 


Same Grid Resources 
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Spectral Volume Grid 


Equivalent Control Volume Grid 


Higher Order Curvilinear Surfaces 





x 

Computing very high order derivatives of metrics is new 





Higher Order Derivatives Improves Spatial Resolution 
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I st Order 
3rd Order 
5th Order 
7th Order 
9th Order 

I I th Order 
1 3th Order 
1 5th Order 
17th Order 
19th Order 
21 st Order 
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Ultra High Fidelity Hermitian Integration, Flux, and 

Interpolation 



Resolves in the “ultra” wavelength range 

- Ultra-short waves (less than two mesh spacings, CKW Tam, 2002) or 
ultra wavelength scaling (10e3 thru 10e-9) 


THE ELECTROMAGNETIC SPECTRUM 
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Advancing in Time 




•Shorthand Notation for Derivatives 

x f _ d a+b+k f(x,y,t ) 

(l . & . fc/ ^ /£* (l f) ^ 

•All time accurate methods effectively do this: 
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Time Derivatives as Space Derivatives — 2 nd Order 
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1 st Order Correction 
2 nd Order Correction 


1 st Order Entropy 
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High Accuracy at Surfaces 
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•Inviscid Surface Boundary Conditions: 


d a (V-n) 

dt a 


d a Q b 

dt a 


Va : (a > 0) 


• Normally, only a=0 is used in CFD 

But to achieve higher order in time at a 
surface must use larger a (Goodrich, 1999) 


Time Cost per Wavelength 
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Efficiency Chart 



Relative Time cost of wave Resolution by Error Requirement 
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Non-equilibrium Model (Gedeon) 
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•Void-Average Velocity Assumption 
may not be good in first 20 layers 
due to density gradients 

•Curvilinear formulation enables 
better surface heat transfer 



Conclusions/Suggestions 
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• Need ultra short wave resolution 

•Use high order parallelizable 


• Easily parallelizable 

methods 


• Multiple Time-Stepping 

•Use DES for oscillating 


• Use High order boundary conditions 

transition modeling 


• High fidelity in space and time 

•Use parallel technology 


• Geometry definition naturally has high 

(Myrinet, Infiniband, Octa- 


fidelity 

Opteron) 


• Look at Spectral element methods, 


Discontinous Galerkin, etc. 


•Extend SAGE/GLIMPS to 3D 

•Extend Manifest to parallel and/or use latest CFD 

•Validate Manifest model with DNS 

•Use regularized regenerator geometry 

•Perform DNS in first 20 layers where non-uniform 
and model the rest 


•Numerical Error 
(Dissipation/Dispersion) 

•Modeling Error- Sliding 
Interfaces, Faceted 
Geometry, Turbulence 
Transition 


